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ABSTRACT 


Transient cavitation produced by pulsed ultrasound in the MHz 
range has been investigated. First, a theoretical background which 
develops bubble motion in general and focuses specifically on transient 
cavitation is presented. Second, the development of an apparatus 
which is designed specifically for studying this type of cavitation is 
discussed. Third, the use of the apparatus to make measurements of 
cavitation thresholds at various frequencies, pulse durations, and pulse 
repetition rates is discussed. Finally, results of the cavitation thresh- 


old measurements are interpreted using theories developed earlier. 
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I. INTRODUCTION 


The term “cavitation” was first used in the late nineteenth cen- 
tury by Sir Charles Parsons to describe the mechanism that caused 
reduced thrust in marine propellers. Lord Rayleigh, in 1917, identi- 
fied cavitation as the cause of erosion in pump components that were sub- 
ject to reduced liquid pressures. Cavitation has come to be defined as 
the rupture of a liquid or a liquid/solid interface due to the reduction 
of liquid pressure [Ref. 1]. Since the theoretical tensile strength of 
pure liquids is exceptionally high, extremely large negative pressures 
(= -1000 bars) should be required to cause the liquid bonds to break. 

Since rupture in real liquids can occur at relatively low negative 
pressures (= -1 bar), it has been postulated that inhomogeneities exist 
within the fluid which cause this reduced tensile strength. Further- 
more, experimental results point to the necessity of these micro- 
scopic inhomogeneities containing gas [Ref. 2]. Thus, cavitation is 
actually the formation of macroscopic bubbles from microscopic bub- 
bles (cavitation nuclei). Acoustic cavitation can therefore be defined as 
the activity of bubbles within a liquid that is subject to an acoustic 
pressure field. 

Microcavitation, or the formation of micron-sized bubbles from 
cavitation nuclei due to the influence of ultrasonic acoustic energy 
(ultrasound), has benefits when used in such applications as ultrasonic 
cleaning and etching, acceleration of certain chemical reactions, and, 


medically, to destroy unwanted cell formations. Microcavitation may, 


however, produce adverse biological effects as well. In particular, the 
improper use of diagnostic or therapeutic ultrasound could possibly 
result in the unintentional destruction of or damage to healthy cells. 
Investigation of microcavitation phenomena and their possibly 
adverse biological effects has produced a significant amount of infor- 
mation on cavitation thresholds. Robert Apfel’s studies of pulsed 
ultrasound cavitation effects [Ref. 3] have summarized the possible 
biological effects in single cells due to excessive sound pressure 
amplitudes. His work has accounted for the effects of frequency, liq- 
uid inertia, and fluid viscosity in the cavitation process. H. G. Flynn 
has done extensive numerical modelling of cavitation phenomena, also 
taking into account viscosity and surface tension as well as compress- 
ibility and heat conduction [Ref. 4]. R. A. Roy, et al. [Ref. 5], and A. A. 
Atchley, et al. [Ref. 6], have produced precision techniques for mea- 
suring transient acoustic cavitation thresholds by detecting the light 
and sound emissions produced during violent bubble collapse. It is 
their research that has become the basis for the research reported in 
this thesis. Although some aspects of cavitation thresholds have been 
previously investigated, there has been little work done to determine 
the time evolution of the cavitation event. The objectives of this 
research are to (1) develop an apparatus in which cavitation can be 
induced and detected simultaneously by two methods— acoustic 
scattering and sonoluminescence; (2) determine transient cavitation 
thresholds; and (3) examine the time evolution of the transient 


cavitation events. 


II. THEORY 


A. BACKGROUND 

The behavior of bubbles under the influence of an acoustic pres- 
sure field is indeed complex. In this study, we will examine the pres- 
sure amplitudes of pulsed ultrasound which will induce the rapid 
growth and collapse of cavitation nuclei. We will vary the frequency of 
the acoustic signal, the number of acoustic cycles per pulse, and the 
rest time between pulses in order to more fully understand the effects 
of each parameter on the pressure amplitude required to cause the 
nuclei to behave in this manner. In addition, we will develop an appa- 
ratus that provides for both acoustic and optical observation of the 
cavitation events. We must therefore be aware not only of the behavior 
of the bubbles when influenced by the sound field but also of their 
propensity for scattering sound and producing visible light. 

In this chapter, we will review the above-mentioned aspects of 
bubble behavior. We will first look at the three regimes of acoustic 
cavitation: stable cavitation, transient cavitation, and gas body activa- 
tion. We will describe the theory that describes bubble activity under 
static equilibrium conditions, under the influence of a continuous 
sound pressure field in which controlled bubble oscillations occur, 
under conditions which result in controlled growth of the bubbles, and 
under conditions which lead to unstable growth and violent collapse. 
A theory will be described that allows one to predict the peak acoustic 


pressure which produces the collapse of a gas bubble. Finally, we will 
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examine the theories which explain both light emission (sonolumines- 
cence) upon bubble collapse and acoustic scattering from multiple 


bubbles. 


B. CAVITATION REGIMES 

It is generally acknowledged that two, and possibly three, cavita- 
tion regimes exist. The focus of our research will be transient cavita- 
tion produced by pulsed ultrasound. As the name indicates, this 
regime is transient in nature; bubbles are likely to exist for less than a 
few acoustic cycles. The results of the measurements to be described 
later indicate that effects of the second regime, stable cavitation, may 
also be important. Stable cavitation occurs when a bubble subjected to 
a sound field oscillates about its equilibrium radius. This regime is 
characterized by bubbles that are relatively permanent; noticeable 
changes in the equilibrium radius occur only after many cycles. The 
third regime (considered by many not to be a true cavitation regime) 
is known as gas body activation, which involves the oscillations of gas 
pockets permanently embedded in cracks, crevices, or within biologi- 
cal systems. As gas body activation does not include actual formation 
and growth of gas bubbles, only stable and transient cavitation will be 


discused further in this chapter. 


C. BUBBLES UNDER QUASI-STATIC EQUILIBRIUM CONDITIONS 
The first stage of any cavitation problem is nucleation, or the for- 
mation of a macroscopically observable bubble from a cavitation 


nucleus. This formation usually involves the loss of mechanical 


11 


stability of the nucleus. Following the developments of Walton and 
Reynolds [Ref. 7], Noltingk and Neppiras [Ref. 8], and Blake [Ref. 9], a 
brief examination will be made of the conditions required for a spher- 
ical bubble to lose its static equilibrium in a liquid. 

Given a bubble of radius Ro in static equilibrium with a liquid of 
infinite extent at a hydrostatic pressure Po, the pressure inside the 


bubble is given by 


P= Po+(R5} (2.1) 


The term (RS is known as the Laplace or Young pressure and is due to 
the surface tension of the liquid, o (the usually minor contribution of 
the vapor pressure has been neglected). Assuming that the liquid 
undergoes an incremental change in hydrostatic pressure to a new 
value P.., the bubble radius will undergo an incremental change to 
some new radius R. Assuming isothermal conditions and an ideal gas, 


we note 


P1V1 = PoVo = constant 


where 
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The assumption of isothermal conditions is justified because the rate 
of evaporation and condensation of vapor at the bubble surface is suffi- 
cient to maintain this condition. After combining these two expres- 


sions and cancelling common terms, we find that 


Pe =(Po+Re|(Fe) = (2) 


Figure 2.1, which is a plot of P. as a function of R, offers insight 
into the stability of bubbles. Assuming there is no gas diffusion into or 
out of a bubble, P.. is the value of the liquid pressure necessary to 
balance the sum of the internal pressure (which tends to make the 
bubble expand) and the Laplace pressure (which makes the bubble 
collapse). P.. is a minimum for a bubble of radius R*. Bubbles of radius 
less than R* are in stable equilibrium, whereas bubbles of radius 
greater than R* are in unstable equilibrium. To see this more clearly, 
consider a bubble of radius R<R*. If R is increased, P.. decreases 
because the decrease in the Laplace pressure | is less than the 
decrease in internal pressure. Consequently, the external liquid 
pressure needed to maintain equilibrium is now less than the actual 
external pressure. As a result, the bubble shrinks back to R. In other 


words, the collapsing pressure dominates over the 
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Graph of Pressure vs. Bubble Radius Showing Stable 
and Unstable Equilibrium States of Bubbles 
Within a Liquid (from [Ref. 7] 
expanding pressure. If the radius had been made smaller, P.. would 
increase corresponding to a dominance of the expanding pressure. 
Again, the bubble returns to R. Therefore, the bubble is in a state of 
stable equilibrium— perturbations in the radius cause the bubble to 
return to its initial radius. By invoking a similar argument for bubbles 
of radius greater than R*, it is seen that the equilibrium is unstable. 
An increase in radius causes P.. to increase due to the dominance of 


the expanding pressure. The bubble would increase in size 


14 


indefinitely. A decrease in radius causes the bubble to shrink to the 
corresponding stable value of R. Because of its role in determining the 
stability of bubbles, R* is called the critical radius. 

Now, consider a bubble of radius Ro < R* in equilibrium at 
pressure Pg. According to Figure 2.1, it is stable. If the pressure is 
lowered to some value Po - Pa, the bubble will grow until it reaches a 
new stable equilibrium radius R' (since R' < R*). If, however, the 
pressure is reduced to Po - Pp, where Pp is the Blake threshold 
pressure, the bubble reaches the critical radius R*. Any further 
reduction in pressure provides no equilibrium value for the bubble 
radius, and the bubble grows uncontrollably. To solve for R*, we must 


minimize P.N. as a function of radius. 


Rewriting equation (2.2) slightly, 
Zio) \ Ses Ya a 20 
P= =(Po +e) RO: [Re] - Fe (2.3) 


then differentiating with respect to R and Setting the result equal to 


zero, we find 


20\,3 |1/2 
3{Po ~ Re RO 


kK = — 9. : (2.4) 


Previously, we considered the Blake Pressure to be a function of 
initial bubble radius. The converse view is to consider the initial 


radius to be a function of Blake pressure. Accordingly, the Blake 
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radius Rg is the initial radius of a bubble that, when subjected to the 
(acoustic) Blake pressure Pp, grows to a radius of R*, becomes 
unstable, and grows uncontrollably. To find the Blake threshold 


pressure, we first rewrite R* (equation 2.4) in terms of Rp: 


20\.3 11/2 
3{Po + Re RB / 


R* = = ae . (2.5) 
Letting 
2023 
K = [Po +R Re (2.6) 


allows R* to be expressed more simply as 


3K\1/2 
R* = 35) : (274 
Or 
K 20 
ers > aa 


Setting Ro = Rp, R= R*, and Po = Po - Pp in equation (2.3), we find 


20 20 





Po = Po -PB=3Re-Ree (Zea) 
or 
4 
Po - Pp = - 3 py. (2.10) 


16 


Rewriting this expression using equations (2.7) and (2.6), we find 


8 36 1/2 
Gmenns = gee —"5g_3_ (oye Wl) 
2(Po + Ra Bs) 
which can be written in normalized form as 
Pp 4 3 XB 1/2 


Z 
where Xp = Robo" 


Summarizing this discussion, a bubble of radius Rg in stable equi- 
librium will lose its stability if subjected to an acoustic pressure of 
amplitude Pg. Upon loss of stability, the bubble will expand to many 
times its initial size; i.e., it will form a macroscopically observable 


bubble. 


D. BUBBLES UNDER THE INFLUENCE OF A SOUND FIELD 

We begin the examination of bubble motion under the influence of 
a sound field with a description of bubble motion by Rayleigh [Ref. 10]. 
He initially considered the collapse of an empty spherical cavity 
located in an infinite mass of incompressible fluid. The hydrostatic 
pressure in the liquid is assumed constant. Rayleigh observed that the 
work done by the collapsing bubble starting from rest had to be equal 
to the kinetic energy of the fluid. Thus, the following holds true: 


Poo = (Ri —~ R3) = =p | v2(r) 4nr2dr (2.13) 


le, 


where r is the distance from the center of the bubble to a point in the 
liquid that is greater than the radius of the bubble R and vi(r) is the 
velocity of the fluid. P.. is the hydrostatic pressure of the liquid and — 
Ro RS) is the decrease in volume of the spherical cavity. The right- 
hand term is the kinetic energy of the liquid, = mv2. Since the fluid is 


incompressible, the change in volume of the liquid is equal to the 


change in volume of the bubble and the following relationship is true: 
4nr2 v(r) = constant. 
R4 
Therefore, 4nr2 v(r) = 4nR? v{R) or v2 (r) = Tz v2(R). Letting v(R) = R, 
it can be shown that 


R4 R2 
ed (2 aa 


Substituting this result into equation (2.13) and performing the inte- 


gration on the right-hand side yields 


P.. = (Ro - R3) = 2n p R3 R2 | (2.15) 


Solving for R2: 


3 
. Ro 
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Since R is equal to it is apparent that dt = = As R decreases, 
however, R is compelled to approach infinity. At this point, Rayleigh 
strives to avoid this obvious flaw by admitting that the initial assump- 
tion of zero internal pressure is invalid. He then assumes that the 
cavity is filled with a perfect gas undergoing an isothermal compres- 
sion. If this assumption is true, Boyle’s law holds and the pressure 
within the bubble will increase until the bubble begins to rebound. If 
fluid viscosity is zero and all other losses ignored, the bubble will con- 
tinue to oscillate between this minimum radius and the initial radius 
Ro. Again invoking the assumption that initial pressure inside the 


bubble is zero or at least some constant. Rayleigh has derived the fol- 


lowing equation governing pressure within the liquid, p, at any point: 


3 3 
4 
p R | Xo R4 | Xo 
Po = a {®_)_a aay 21.7) 


The work of Plesset, Poritsky, and Noltingk and Neppiras con- 
tributed to the final formulation of an equation that describes the 
motion of a bubble under the influence of an acoustic field. Plesset 
(Ref. 11] extended Rayleigh’s analysis, as did Poritsky [Ref. 12] and 
Noltingk and Neppiras (Ref. 13]. In 1976, Lauterborn [Ref. 14] com- 
bined these results with previously described results to form the final 


characteristic equation governing bubble motion, the RPNNP equation. 
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Lauterborn so named this equation in honor of the pioneering 
work done by Rayleigh, Plesset, Noltingk, Neppiras, and Poritsky. In 
review, he set the work done by the liquid at infinity minus the work 
done by the liquid layer immediately adjacent to the bubble equal to 
the kinetic energy of the liquid surrounding the bubble: 


f (PL - P..) 4nr2 dr = 2npR3R2. (2.18) 
Ry 


where Py is the pressure in the liquid layer immediately surrounding 
the bubble. The reader should note that P. in this development is 
varying with time: 


where Po is a static pressure about which P(t) varies. Therefore, from 


Reference 14, Py, can be expressed as 


pL=(Po+Re-Pv\(qe) -2R- UR (2.19) 


where x is the polytropic exponent of the gas inside the bubble, up is 
the shear viscosity of the liquid, and Py is its vapor pressure. 
Differentiating (2.18) and substituting (2.19) for PL and Px = Po + Pit), 
we find that 


RR + >R2 = 1\(p +e-Pvl ey 20 SUR Pl) eo. 


The following assumptions are implicit in this comprehensive 

equation [Ref. 7]: 

1. The bubble remains spherical at all times; 
. Uniform conditions exist within the bubble: 
. The acoustic wavelength, 4, is much greater than R; 
. No body forces are present; 
Bulk viscosity is ignored; 
. The liquid is incompressible; 


. Gas content of the bubble is constant— no diffusion occurs: 
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. Vapor pressure, Py, is constant during the bubble’s motion. 

Although many simplifying assumptions have been made in deriv- 
ing equation (2.20), it is a better-than-reasonable description of the 
dynamics of a bubble subjected to an ultrasonic sound field [Ref. 7]. In 
addition, the RPNNP, or simply the Rayleigh-Plesset, equation is the 


foundation of the theories discussed in the next two sections. 


E. STABLE BUBBLE GROWTH DUE TO AN ACOUSTIC PRESSURE 
FIELD 


Stable cavitation and rectified diffusion, although neglected in the 
RPNNP equation, are two related and interesting processes which may 
be important in the time evolution of transient cavitation events. The 
stable growth of cavitation nuclei to sizes which allow transient growth 
is a likely mechanism for producing transient cavitation. Rather than 


undertake a rigorous mathematical development of rectified diffusion 
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as in Reference 15, a qualitative description of this phenomenon will 
be presented here. 

The importance of rectified diffusion, or incremental enlarge- 
ment, was first noted by Harvey, et al., in 1944 [Ref. 16]. Blake [Ref. 9] 
was the first to attempt the derivation of a mathematical model in 
1949. Two effects, the area effect and the shell effect, appear to work 
in concert to produce bubble growth. 

1. The Area Effect 

When a bubble contracts, the concentration of gas within the 
bubble increases. Since gas diffuses from regions of high concentra- 
tion to those of low concentration, the gas diffuses out of the bubble. 
When the bubble expands, the concentration of gas decreases, and the 
concentration gradient reverses, the result being a diffusion of gas into 
the bubble. The diffusion rate is proportional to the area through 
which gas may flow. Therefore, more gas enters the bubble when it is 
expanded than leaves when it is contracted. The net result is an influx 
of gas during one complete oscillation of the bubble. 

2. The Shell Effect 

The rate of diffusion, as noted above, is proportional to the 
concentration gradient. Crum [Ref. 15] considers a spherical shell of 
liquid surrounding a bubble. He argues that, when the bubble contracts, 
the shell expands, with a consequent reduction in gas concentration near 
the bubble surface. This increases the concentration gradient in the 
outward direction over that discussed above, so the diffusion of gas out 


of the bubble is even greater. Conversely, when the bubble expands, 


De 


the shell volume decreases, the gas concentration increases, and the 
concentration gradient into the bubble increases. The overall result is 
the enhancement of the area effect. 
3. Threshol ion 
There exist both high- and low-frequency thresholds for 
rectified diffusion. Crum has formulated expressions for both. Of 


interest here is the high-frequency approximation [Ref. 15]: 


2 2 | al 20 eH 
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DPM) 
Ci 26 ( 
(3+4K) Cy —-K/ 1+ Ra Pa) 
where 

Pa = pressure threshold for rectified diffusion 

K = 

2 
B2 = p w2Ro/3 P.. 
b = 4oy/3 P. (Damping term) 
2 ] 40 
oon = (—=—7 (s P.. +f) (small amplitude resonance frequency 
PRo Squared) 


Co = equilibrium concentration of gas in the liquid 
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C; = concentration of dissolved gas in the liquid far from the 


bubble 
1 = viscosity of the liquid‘ 
® = angular frequency of the acoustic field 


Ro = equilibrium bubble radius 


po = liquid density 
P,. = ambient pressure 
o = surface tension of the liquid 


F. BUBBLES DRIVEN TO UNSTABLE GROWTH AND VIOLENT 
COLLAPSE 


The final step in development of the theory which defines bubble 
behavior in a sound field is taken from research by R. A. Apfel and H. G. 
Flynn. In Reference 3, Apfel indicates that the possibly adverse bio- 
logical effects of acoustics are linked to the growth and collapse of 
submicron cavitation nuclei. A nucleus “that grows to 5 um in one 
acoustic cycle will, upon collapse, deposit in the size scale of one bio- 
logical cell approximately 300 million electron volts (MeV)” of energy 
[Ref. 3]. This energy, in the form of heat, results in internal gas tem- 
peratures in the range of 1000-5000 K. At such high temperatures, 
effects such as free radical formation and recombination can occur. 
These effects can result in the production of visible light that can be 
detected optically (this will be discussed in detail later in this 
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chapter). Apfel’s work culminates in a formula that predicts the peak 
acoustic pressure which produces a cavitation event (collapse) given 
the initial bubble radius Ro and maximum internal gas temperature 
9000 K with acoustic frequency, density, viscosity, and surface tension 
of the liquid as parameters. He develops the relationship between the 
maximum collapse temperature T’ and normalized acoustic pressure p 


I 
2 Eo (Pa is peak acoustic pressure; Po is ambient pressure) by relating 


p to Rmax (maximum bubble radius) and Rmax to T’. Although a rigor- 
ous development of these two relations will not be conducted here, 
Figure 2.2 shows readily how bubble size corresponds to applied 
acoustic pressure. The reader is referred to Reference 2 for a more 
rigorous development of T’ vs. Rmax and Rmax vs. p. It is necessary, 
however, to examine Figure 2.2 carefully to determine the character of 
the Rmax vs. p relationship. The maximum bubble radius can be 
thought of as being composed of two parts: Rj, which occurs when 
AP, the pressure difference across the bubble, is negative (outward 
force); and Rmax, which occurs when AP is positive, yet the bubble 
continues to expand due to outward liquid momentum. 

The time over which maximum bubble radius is computed is also 
noted in Figure 2.2. The time t, during which AP is negative, is the 


algebraic sum of tg, ti, and At. 


t=t2-t, -At, (2.22) 
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Figure 2.2 


Schematic of Bubble Radius as a Function of 
Time When Subjected to a Single Acoustic Cycle (from [Ref. 3]) 


where 
} = 
to mae a, (2.23) 
ts - 218) (2.24) 
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and 


2 
at = 2 Ro Vf o5+ 48. (2.25) 


Ee 
Padenen | — Bo and pp is the normalized Blake threshold pressure given 


by equation (2.12). The value of t is therefore 


em , /2te-0) ny Ce 2 = pee Of 
| p p mito 2 ba ip 1) Po (p-1) 
= 7 


(2,26) 


The maximum bubble radius can now be estimated by 


- 1/3 
OSE Re a Cue PHT |[1 +2 - | (2.27) 


PR3 
From the perfect gas law, =;— = constant, and its adiabatic form PR*’ 


= constant, it can be shown that 
Peet (y= 1) Ta (2.28) 


From equations (2.26) and (2.27), above, we find 


T =To v=o +2R5 EP [13 to (2.29) 


It is now convenient to find frequency (embedded in 7), resulting in: 
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i (2.30) 
S20 _ au. /_ 2 9.46 
‘a Ms Ro Vp Pox 
3+ ] 
where 


X=Pp+p-2+V(p-1) (p - pp). 


In Table 2.1 from Reference 3 and Figure 2.3 from Refer- 
ence 4, the effects of initial bubble radius and frequency in transient 
cavitation threshold can be seen. It is important to note that fre- 
quency and internal temperature have little effect on the threshold 
pressures of smaller bubbles, but become increasingly important as 
initial bubble radius gets larger. This approximate formulation is not 
meant to be a substitute for more exact numerical calculations such as 
are found in Reference 4, but it does give a reasonable idea of the 
manner in which the parameters of frequency, viscosity, and initial 


bubble radius affect the transient cavitation threshold of a liquid. 


G. SONOLUMINESCENCE 

Since one of the goals of this thesis is to develop both optical and 
acoustical cavitation detection schemes, it is important to present 
some theories concerning sonoluminescence. Following the review of 
Walton and Reynolds [Ref. 7], two primary theories as to the origin of 


sonoluminescence have been advanced. 
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TAB TE 
PRESSURE THRESHOLDS PREDICTED BY EQUATION (2.30) 


Pressure (bar) Pressure (bar) 


T° K. f(MHz) Ro = 0.l um Ro = 0.5 um 
2500 1.0 6.4 Zee 

9000 6.4 2.1 

2500 Zo 6.4 2.8 

5000 6.4 oo 
2500 5.0 6.4 4.2 

9000 6.4 Daal 


Pulse wiath 10°78, 


* Pressure threshoid, Py (bars) 





initial Radius, Rp_ (em) 


Figure 2.3 


Threshold Pressure P; (in Bars) for Transient Cavitation 
as a Function of Initial Radius Ro Over a Range of 
Frequency fa, at a Pulse Width of lus [Ref. 4] 


pas, 


Noltingk and Neppiras [Ref. 13] proposed that sonoluminescence 
is the result of black-body radiation from the gas contained within the 
bubble, which theoretically rises to extreme temperatures as the bub- 
ble undergoes a rapid adiabatic compression during its collapse. 

Griffing [Ref. 17] proposed that the formation of highly reactive 
free radials in the liquid during bubble collapse were responsible for 
rapid chemical changes in the fluid, resulting in the emission of light. 

For the most recent examinations of sonoluminescence, the 
interested reader is referred to References 18 through 22. The pre- 
ponderance of evidence points to the formation of free radicals that 
subsequently play a role in sonoluminescence. These free radicals, or 
“reactive molecules (or atoms) with unpaired electrons” [Ref. 23] may 
be formed as a result of the dissociation of molecules subjected to the 
high temperatures produced by a rapid adiabatic compression. Often, 
if a free radical (with its unpaired electron) reacts with a molecule 
having its full complement of paired electrons, the result is another 
free radical. This chain of events will be self-sustaining until all of the 
Suitable material has been reacted with, the free radicals have com- 
bined with each other to form stable molecules, or all of the free 
radicals have collided with the wall of the container, where they are 
removed from the process. According to Griffing, the emission of 
light may be due to one of these free radical reactions that happens to 
be of a chemiluminescent nature. This hypothesis has been confirmed 


by a number of subsequent experiments [Refs. 18-22]. 
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H. ACOUSTIC SCATTERING FROM BUBBLES 

The other detection method that we hope to use to investigate 
transient cavitation is the scattering of an acoustic signal (either the 
source signal or an alternate lower-frequency active signal) from 
bubbles generated by such events. This technique is expected to take 
advantage of the exceptionally large scattering cross-section of 
resonant bubbles. Since our measurements do not require numerical 
solution of scattering parameters, however, a qualitative explanation 
rather than a rigorous mathematical development will be given here. 
This theory is taken from Clay and Medwin [Ref. 24]. 

A small, rigid sphere whose density and elasticity are greater than 
that of the medium backscatters a very small amount of the acoustic 
energy when ensonified by acoustic frequencies where the wavelength 
is much greater than the sphere’s radius. The acoustical cross- 
section of the sphere is, in fact, much smaller than the geometrical 
cross-section because of diffraction of the incident sound field around 
the small, rigid body. 

Bubbles, however, differ from rigid spheres in two ways: (1) the 
specific acoustic impedance (pc) of the gas within the bubble is much 
less than that of the water, and (2) bubbles resonate. Bubbles are 
extremely effective absorbers and scatterers of sound at frequencies 
near resonance frequencies. Scattering cross-sections at resonance 
are on the order of 400 times their geometrical cross-sections. (See 


Figure 2.4.) 


oul 


2 
0, /Na 


Bubbie 


Rigid 
spnere 





Figure 2.4 
Ratio of Scattering Cross-Section to Geometrical 
Cross-Section for a Bubble and a Rigid Sphere (from [Ref. 24)]) 
Transient cavitation events studied in this investigation will result 
in the generation of many bubbles of varying sizes. The effects of 
multiple scattering must be considered. When scatterers are widely 
spaced, the scattering cross-section per unit volume is simply the 


sum of the individual scattering cross-sections of the bubbles within 
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that volume. The bubble “cloud” expected as a result of the transient 
cavitation event will be much more densely packed, however, which 
will cause the individual scattered fields of the oscillating bubbles to 
interact. The result is that the resonance curve of a single bubble 
shown in Figure 2.4 is broadened and the Scattering cross-section of 
an array of bubbles is less than the sum of the individual cross- 
sections. A cloud of bubbles— bubbles packed randomly in a three- 
dimensional array— acts as a pressure release surface at the face of the 
cluster, thus reflecting acoustic energy at this interface. It is this 


effect that is expected to reveal the presence of acoustic cavitation. 
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Ill. EXPERIMENTAL APPARATUS 


A. BACKGROUND 

The apparatus with which these cavitation measurements will be 
made is based upon a system developed by A. A. Atchley, L. A. Frizzell, 
and R. E. Apfel [Ref. 6]. Our apparatus differs from theirs in three 
major areas: (1) cavitation detection can be accomplished by two 
methods: optically (via sonoluminescence) and acoustically, either 
individually or simultaneously; (2) this apparatus will rely on random 
bubble nuclei attaining threshold radii vice providing an artificially 
introduced uniform cavitation nucleus to the sound field; and (3) this 
apparatus is computer-controlled. The system is made up of two 
major subsystems: the fluid management system and the electronics 
system. The fluid management system will distill, filter, and degas 
ordinary tap water for use in the cavitation enclosure. The electronics 
system will provide the capability to induce, detect, and record tran- 


sient cavitation events, all under computer control. 


B. CAVITATION ENCLOSURE AND LIGHT-TIGHT ENCLOSURE 
The cavitation enclosure (Figure 3.1) is a 25 cm by 15 cm by 15 
cm box made of 1.2 cm thick plexiglass. Plexiglass was chosen as the 
primary structural material because it is sufficiently strong to support 
installed transducers and fluid inlet and outlet connections, yet it 
allows the optical detection apparatus to “see” what occurs in the 


insonified region of the system. One end of the chamber holds a 
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Figine so: | 


Cavitation Enclosure 





Figtire 37 


Light-Tight Enclosure and Photomultiplier Tube Assembly 
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Harisonic Labs focused transducer that supplies the cavitation genera- 
tion signal. A wedge and two vertical walls (also plexiglass) have been 
installed at the opposite end of the enclosure to scatter the incident 
sound and thus minimize reflections that would interfere at the focal 
region. The source transducers (rated at 1.0 MHz and 2.25 MHz) have 
focal lengths of 7.6 cm. Two additional walls with holes cut to allow 
passage of the focused sound beam are installed on either side of the 
focal region, again to minimize interference from reflections. The 
mounting for the source transducer allows for rapid changing of ele- 
ments so that the other source frequency can be examined. Two 
Panametrics V301 transducers are permanently installed in the bot- 
tom and one side of the enclosure in such a manner that their beam 
pattern intersects the focal region of the source transducer. The top 
of the enclosure is removable to allow access to the interior. The 
cavitation enclosure is placed in a light-tight box that has connections 
for fluid inlet and outlet as well as BNC connections for source and 
receive transducer signals (Figure 3.2). The photomultiplier tube 
(PMT) assembly is mounted on one side of the light-tight enclosure so 
that it also can “see” the focal region for sonoluminescent indications 
of cavitation. The photomultiplier tube’s sensitive circuitry is pro- 
tected by a mechanical shutter that can be lowered to block light and 
by an interlock switch in the lid of the enclosure that interrupts the 


high-voltage power supply to the PMT when the lid is removed. 
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C. FLUID MANAGEMENT SYSTEM 

The fluid management system (Figure 3.3) consists of a four-ele- 
ment deionization/filtration system, a bag-type 0.2 um filter, a small 
variable-speed electric centrifugal pump, an ultraviolet sterilizer, and 
an Erlenmeyer flask with connections for system evacuation and argon 
infusion. The Erlenmeyer flask is also used as a fluid reservoir and 
expansion tank. The filtration/deionization system produces the equiv- 
alent of quintuply distilled water and removes all ionized minerals, 
silicon, and free carbon dioxide down to four parts per billion. The 
bag filter ensures that no particles larger than 0.2 um are introduced 
into the system. These two elements are used only when the system 
is being filled. When the system is fully charged, the water is circu- 
lated through the cavitation enclosure via a closed loop consisting of 
the pump, the U/V sterilizer, the enclosure, and the expansion reser- 
voir. Argon is bubbled into the fluid to help enhance the sonolumines- 
cent reaction. The system is simultaneously evacuated (15 inches Hg 
vacuum) in order to remove dissolved gases. The evacuation is contin- 
ued until visible air bubbles are removed (usually about one hour). The 
circulation system is then secured and the sample is returned to 
atmospheric pressure, where it remains quiescent for one-half hour 
before data accumulation begins. After about an hour of data accumu- 
lation (usually 20 to 40 data points), the system must be drained, refilled, 
and reprocessed. Refilling the system is necessary because the fluid 
becomes “tired” after multiple cavitation events, causing higher-than- 


normal thresholds (possibly due to fewer remaining cavitation nuclei). 
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D. CAVITATION GENERATION AND DETECTION SYSTEM 

The electronics system consists of the cavitation generation sys- 
tem and the cavitation detection system (Figure 3.4). The heart of the 
cavitation generation system is the Hewlett-Packard 9000 Series 300 
personal computer, which controls the output frequency, pulse width 
(duty cycle), pulse repetition rate, and acoustic pressure amplitude of 
the source transducer via the Hewlett-Packard 3314A function gener- 
ator. The output of the function generator is boosted by an Amplifier 
Research 100-A15 Power Amplifier and is then applied to the source 
transducer. A printout of the control program is contained in 
Appendix A. 

The cavitation detection system consists of two subsystems: the 
acoustic detection system and the photomultiplier tube (PMT) assem- 
bly. The acoustic detection system is made up of the two Panametrics 
V301 transducers installed in the cavitation enclosure and the associ- 
ated elements for either passive or active (via scattering of a cw signal) 
detection of cavitation events. In the passive mode, the output of one 
of the detection transducers is boosted by a small RF amplifier and is 
displayed on the Tektronix 2245A oscilloscope. The trace of this 
signal on the oscilloscope consists primarily of a bundle of acoustic 
pulses that have either been diffracted or scattered to the detection 
transducer face. The bundle is delayed from the source pulse by 
approximately 100 microseconds, which corresponds to the time the 
pulse would take to make the 15 cm journey from the source to the 


focal region and then to the detection transducer face. The 
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system operator must monitor the signal in the vicinity of this time 
delay for a transient that would indicate an increase in the scattered 
source signal, thus signifying the occurrence of a cavitation event. The 
operator must then manually interrupt the computer ramp function 
for the computer to calculate and display the threshold pressure. 

In the active mode, the second Panametrics transducer is driven 
at a continuous frequency in the 100 kHz-500 kHz range. The result- 
ing active signal is received by the detection transducer along with the 
source signal. The output of the detection transducer is passed 
through a low pass filter to remove the cavitation generation signal 
while passing the active detection signal. A sharp increase in the 
amplitude of the received signal indicates scattering from the bubble 
cloud that is formed by a cavitation event. The photomultiplier 
assembly and associated circuitry provide optical detection of the 
sonoluminescence generated by a transient cavitation event. The 
pulse from the PMT is pre-amplified and sent to an ORTEC 9302 
pulse height discriminator, which determines the sensitivity of the 
mechanism. The discriminated signal is then sent through a Canberra 
2071A counter. Because of background light noise within the light- 
tight enclosure, there is always some pulsed output from the PMT. 
Since transient events nearly always occur in groups, events will likely 
produce numbers of pulses significantly higher than background. At 
the time of this writing, attempts to detect cavitation via this method 


have proved unsuccessful. 
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E. TRANSDUCER CALIBRATION APPARATUS 

In addition to the apparatus used for the generation and detection 
of cavitation, the apparatus used for the source transducer calibration 
is sufficiently involved to warrant further description. Calibration 
involves measuring the transverse and axial beam pattern of the source 
transducer and determining the transducer’s total acoustic power 
output by force target deflection. 

The axial and transverse beam patterns of both the source and 
detection transducers were found using a small PVDF probe 
hydrophone designed and built specifically for this use. The 
construction and calibration of this hydrophone were accomplished to 
fulfill requirements for a final project in a transducer theory and 
design course. Pertinent portions of the report submitted on the 
project are contained on Appendix B. The PVDF hydrophone was 
mounted in a micrometer jig that allowed precise measurement of its 
position relative to the transducer face (Figure 3.5). 

The force target deflection measurement was used to determine 
peak output frequencies and total acoustic power output of the source 
transducers (Figure 3.6). A description of the theory behind the force 
target procedure follows in Chapter IV. In this procedure, the source 
transducer is driven by the function generator via the power amplifier 
at the same gain setting as will be used in the threshold experiment. 
The target is a small plexiglass plate, suspended at a 45-degree angle 
by four pieces of two-pound test monofilament fishing line of equal 


length (approximately .04 mm diameter, 225 wgrams/cm linear 
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PVDF Hydrophone Determination of Source Transducer Beam 
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Measurements 
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density). The plexiglass target is rectangular in shape, of dimensions 
3.8 x 5.1 x 0.7 cm. It has been carefully weighted using lead strips to 
achieve a weight of 1 g in water (14.7 g in air). A cavity approximately 
O.3 cm deep is machined into the target and covered by a clear, thin, 
stiff plastic to form a pressure-release surface. This ensures that the 
radiation force is almost entirely transformed into a deflecting force. 
Essentially no acoustic energy is reflected back to the driving 
transducer, but rather toward the bottom of the test tank. The 
deflection was measured by a graduated telescope mounted to the 


micrometer jig described above. 


Ad 


IV. EXPERIMENTAL PROCEDURE 


A. BACKGROUND 

Prior to data taking, it is necessary to determine the performance 
characteristics of the transducers involved. This involves reciprocity 
calibration of the two Panametrics V301 transducers used for the 
acoustic detection of cavitation events. Also, once the source 
transducers’ true peak output frequencies are determined, it is 
necessary to determine their axial and transverse beam patterns. 
These transducers are used to generate the ultrasonic field that 
induces cavitation. The axial beam patterns are used to verify the 
transducers’ focal lengths while the transverse beam patterns are 
integrated to yield the effective area, Aeff, which relates peak intensity 
and total transducer output power. 

Having determined the beam patterns, it is necessary to perform 
a radiation force-target deflection measurement of the Harisonic 
transducers in order to calculate transducer output power. Having 
experimentally determined Aeff and total power, the output peak 
intensity and consequently the output pressure characteristics can be 
accurately determined. 

Upon completion of transducer source strength calibration, cavi- 
tation samples are prepared and threshold measurements are made in 


order to establish a data base. 
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B. RECIPROCITY CALIBRATION 

The technique employed to find the sensitivities of the Panamet- 
rics V301 transducers is identical to that described in Reference 25. 
See Figure 4.1. The distance between transducers is 20.2 cm. This 
distance is well into the far field of the transmitter, which is found 


using the following approximation: 





in zz (4.1) 
where 
r = acoustic wavelength 1.5 mm (assuming c = 1481 m/s for 
fresh water and f = 986 kHz), 
L = diameter of V301 transducer = 3.2 cm, 
I'min = Minimum distance to the acoustic far field (in cm). 


By employing the appropriate equations established in Reference 25, 
graphs of the transducers’ sensitivities can be constructed. See Figure 


4.2. The V301 transducers with serial numbers 596 and 598 were used. 


C. AXIAL BEAM PATTERN MEASUREMENT OF SOURCE 
TRANSDUCERS 


Using the polyvinylidene fluoride (PVDF) hydrophone discussed 
above (see Appendix B) and a precision measuring jig (see Figure 3.5), 
the axial beam patterns of the two Harisonic source transducers are 


determined. A CW sinusoid at the transducer’s rated operating 
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(a) 


(b) 


Figure 4.1 
Reciprocity Calibration of Transducers (from [Ref. 25]) 
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Figure 4.2 


Sensitivity Calibration Results of 3 Panametrics V301 Transducers 
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frequency is generated by the HP 3314A function generator and 
passed through the Amplifier Research model 100A15 power 
amplifier and finally to the test transducer. The amplifier gain is set at 
200 for a constant 15 VAC (peak-to-peak) output. The PVDF 
hydrophone is inserted into the mounting bracket attached to the 
measuring jig. The active element of the hydrophone is oriented such 
that the hydrophone face moves along the transducer’s longitudinal 
axis as the micrometer knob is rotated. By recording the received 
voltage from the hydrophone and the position from the measuring jig’s 
micrometer Knob, an axial beam pattern can be determined and the 
transducer focal length verified. The Harisonic transducers used in 
this series of experiments were found to have focal lengths of 7.29 and 
7.37 cm. There is not a large enough difference between these values 
and the rated focal length of 7.62 cm (3 in.) to warrant redesign of the 
cavitation tank or relocation of the transducers used for acoustic 


detection of cavitation. 


D. TRANSVERSE BEAM PATTERN MEASUREMENT OF SOURCE 


TRANSDUCERS 
1. Theory 


It is necessary to relate peak intensity and total acoustic 
power in order to establish a relationship between transducer input 
voltage and transducer output pressure. In order to establish this 
relationship , we define an effective beam area Aeff by 


PTOT 
IPEAK ’ 





Aeff = (4.2) 
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where IpgfaKk = peak intensity (W/m?) and Prot = total output power 
(W). Since I a p2 a V2, I = K'V2 where K' is a constant: p is acoustic 
pressure and V is the corresponding hydrophone ouput voltage. This 
voltage, a function of distance from the center of the transducer face, 
is determined by a transverse scan of the focal region. At the focal 


2 
point, Ippak = K’'Vpgax. Substituting this expression into equation 


(4.2) gives 
ETO TE 
Aeff = a a (4.3) 
K'VPrAK 


It is known that 


Prot =! IdS; dS = incremental area normal to the beam. 


The integral is taken over the entire cross-sectional area of the beam 


and can be approximated (in polar coordinates) by 


716 
PTOT = | | K'V2rdrd@ = 2nK' | V2rdr (4.4) 


where r = radial distance from transducer axis and a is some radial 
distance within which is contained practically all of the beam energy. 
In practice, a is chosen equal to the radius of the transducer face. 


Therefore, 


ank'| V2rdr an] V2rdr 


Aete = ——3——_ = — —— (4.5) 
RGU RAK VPEAK 
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Finally, since Agff can be determined from the transverse beam profile, 
it is necessary only to calculate Pror (as described in a subsequent 
section) in order to calculate the peak intensity since IpgaK = 
Prot /Aeff- 

Further, since peak pressure, Pprax, iS related to peak 
intensity by the following relationship (assuming either plane or 


spherical waves), 
PPEAK = (poclpgeax) !/2, (4.6) 


PPEAK Can be rewritten in terms of Aeff and Prot: 


oc 
ppeaK = (Aon ot! (4.7) 
2. Procedure 


The procedure used for measuring transverse beam patterns 
is very similar to that described in the previous section. The jig is 
positioned in a manner such that the hydrophone travels in a line 
perpendicular to the longitudinal axis at a distance from the 
transducer face equal to the focal length. Accurate position measure- 
ments (+10 um) can be made by recording the micrometer knob 
readings. Concurrently, hydrophone output voltages are recorded 
from the oscilloscope display, resulting in the transverse beam pat- 
terns of Figures 4.3 and 4.4. Plots of V2r vs. r are included as Figures 


4.5 and 4.6. Only the right half of each beam pattern was used to 
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Figure 4.3 


Transverse Beam Pattern of 1.0 MHz Transducer 
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Figure 4.4 


Transverse Beam Pattern of 2.25 MHz Source Transducer 
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Figure 4.5 


V2r vs. r (Integrand of Equation (4.4)); 
1.0 MHz Source Transducer 
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Figure 4.6 


V2r vs.r (Integrand of Equation (4.4)); 
2.25 MHz Source Transducer 
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be symmetrical. This is not exactly true, but no appreciable error is 
introduced by making this assumption. Integrations of these plots are 
performed by using a Macintosh Plus personal computer and DBGPLOT 


software. The results are plotted in Figures 4.7 and 4.8. 


E. DETERMINATION OF TRANSDUCER PEAK OUTPUT FREQUENCY 

In order to optimize the output from the source transducers, it is 
necessary to know their peak output frequency. Using the deflection 
apparatus illustrated in Figure 3.6, the rectangular target is suspended 
in the focal region of a source transducer. A sinusoidal signal at the 
rated operating frequency is applied to the transducer via the power 
amplifier. A constant input voltage of 15 volts (peak to peak) is main- 
tained (monitored on the oscilloscope) while frequency is varied above 
and below the manufacturer’s quoted peak output values. Target 
deflection is measured for the various test input frequencies by noting 
the change in target position with a graduated telescope mounted on 
the measuring jig used previously. Plots of target deflection versus 
transducer frequency are reproduced as Figures 4.9 and 4.10. The 
transducers specified as 1 MHz and 2.25 MHz peak output frequencies 
were found to have true peak output frequencies of 986 kHz and 2.17 
MHz, respectively. 
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Figure 4.7 


Integral of V@rdr for 1.0 MHz Transducer 
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Figure 4.8 
Integral of V2rdr for 2.25 MHz Transducer 
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Figure 4.9 


Target Deflection versus Frequency for Determination 
of Peak Operating Frequency of 1.0 MHz Source Transducer 
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Figure 4.10 


Target Deflection versus Frequency for Determination 
of Peak Operating Frequency of 2.25 MHz Source Transducer 
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F. RADIATION FORCE-TARGET DEFLECTION MEASUREMENT 
1. Theory 

The method used to determine the peak output pressure of 
the Harisonic transducers is a variation of the method described by 
Dunn, et al. [Ref. 26]. Rather than using a sphere, as in Dunn’s work, 
the rectangular plexiglass target is employed. 

A rectangular target is used rather than a sphere because, as a 
sphere deflects, it tends to move out of the sound field, while a rela- 
tively large plate is less susceptible to this effect. In addition, there is 
no need to calculate the form function. 

Referring to Figure 4.11 (reproduced from Reference 26 with 
some clarifying remarks), it is a relatively simple matter to construct 
the appropriate free-body diagram and the equations governing the 


static equilibrium: 
Fx = -Fr + Tcos@ = 0, 
LFy = -mg + Tsin@ = 0. 


The above equations imply that 


Fy mg 
ta cos8@ ~~” _sin@’ (4.8) 
f Eo H ‘arQus = If L>>h, tan@ = L/d 
Therefore, Fr = and" owever, tand = —j— >>h, tané = , 
Therefore, 
_ mg _mg_ mgd 
d 
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Figure 4.11 


Free Body Diagram of Force Target Deflection 


From References 25 and 27, it can be shown that the excess pressure 


De = pc2s, where p = density of the medium, s = condensation, and c = 
phase speed of the acoustic waves. For a high-intensity sound field, 


variations in density are significant, so p=po(1 +s). By substitution, 
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De = pc*sS = po(l] + s)c2s = poc*s + poc2s2. If s = socosat, the time- 


averaged excess pressure, <De>t, can be rewritten as follows: 


2 
<De>t = Poc2So <cosmt> + poc2So <cos2wt>. 


The first term on the right-hand side reduces to zero, but the second 
term (or radiation pressure) can be expanded by employing the 


appropriate trigonometric identity. 


Z 2 1 
Poc2Sp <cos2e::> = poc2So - 5 [1 + <cos2@t>] 


Zz Pd 
Poc*Sp = poc*So 
=— 3 + Zz <cos2at>. 








Z 
, poc2sg 
This reduces finally to poc?so <cos*mt> = —y—— . However, De = = 





where I = intensity. Recalling that power, Pror =!1- A where A = area, 
it can be seen that the relationship Prot = <DPe>tcA is also true. 
Fy 


Setting <pe>t:- A=Fy,;, Prot = A cA = F;- c where c = 1481 m/s since 


these tests are conducted in fresh water. By substitution, 


mgc 
Prot ===: d (4.10) 
where oo is a constant and d is the measured target deflection. 


By substituting the appropriate values for m, g, c, and L, 
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Daan = (0:99 x 10-3 kg) (9.8 m/s?) (1481 m/s) e 
TOT = 14.05 x 10-2 m oo 


= (1.02 x 102 ot d: d in meters. 


This equation is used in the radiation force-target deflection mea- 
surements. Upon completion of such measurements, all ingredients 


necessary to characterize transducer output pressure will have been 


determined. 
2. Procedure 


A CW sinusoid (at the transducer’s peak output frequency) 
generated by the HP 3314A function generator is sent through the 
Amplifier Research model 100A15 power amplifier and finally to the 
Harisonic transducer undergoing test. The deflection of the plexiglass 
target is measured as the transducer input voltage (as monitored on 
the Tektronix 2445A oscilloscope) is varied. Plots of target deflection 
vs. input voltage squared are included as Figures 4.12 and 4.13. The 
slopes of these plots are proportional to output power for a given input 
voltage. Note that the “mirror” images of the recorded data points are 
plotted as well. This ensures that a line determined by linear 
regression will pass through the origin, resulting in zero power output 
for zero input voltage, as is dictated by common sense. 

Having determined all of the terms necessary to characterize 
output pressure, the pressure as a function of input voltage can be 
plotted. See Figures 4.14 and 4.15. Again, the “mirror” images are 


plotted as well in order to drive the linear regression through the 
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Figure 4.12 


Target Deflection vs. V2 for 
Source Transducer Calibration 1.0 MHz 
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Figure 4.13 


Target Deflection vs. V2 for 
Source Transducer Calibration 2.25 MHz 


60 


Transducer Output Pressure (bar) 


Transducer Output Pressure (bar) 
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Figure 4.14 


Output Pressure Amplitude vs. Input Voltage 
for 1.0 MHz Source Transducer 
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Figure 4.15 


Output Pressure Amplitude vs. Input Voltage 
for 2.25 MHz Source Transducer 
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origin. Note that, in Figure 4.15, the relatively large deviations from 
the linear regression at lower voltages are of no concern since the 
transducer is not operated in this region. The relationships between 
the peak-to-peak transducer input voltage and the peak output 
pressure are included below. 


For the 1 MHz transducer: 


PPEAK (bar) = 0.0522 Vin (mV). 


For the 2.25 MHz transducer: 


PPEAK (bar) = 0.0403 Vn (mV). 


G. DETERMINATION OF TRANSIENT CAVITATION THRESHOLDS 

All threshold data were taken using the passive acoustic (scatter- 
ing of the source signal from the bubble “cloud” produced by a tran- 
sient event) detection scheme. This method was the first to show 
acceptable reliability. As the data required was substantial and excep- 
tionally time-consuming to take, we concentrated on data accumula- 
tion via this method rather than sacrifice more time refining the other 
two detection methods. The HP 9000 series 300 computer was pro- 
grammed to allow user input of desired frequency, pulse repetition 
frequency (PRF), pulse duration, and pressure amplitude sweep para- 
meters (via input of function generator starting amplitude, voltage 
amplitude increase per step, and total number of amplitude steps). 


Data was obtained at 986 kHz and 2.17 MHz using 200 Hz, 500 Hz, 
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ec OOCszee kh Pilse duranonssolela, 2, 510, 15, 20, 25, 30, 35, 
and 40 usec were examined. All voltage ramps were started at 300- 
400 mV peak input voltage to the power amplifier, which results in an 
initial peak output pressure from the source transducers of approxi- 
mately 3-4 bar. Incrementally increasing pressure steps of .1 bar 
were held for 1.5 seconds. The ramps were continued until a cavita- 
tion event was detected. Twenty threshold measurements were taken 


at each pulse duration, PRF, and frequency combination. 
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A. RESULTS AND DISCUSSION 

Figures 5.1 and 5.2 are graphs of threshold pressure versus pulse 
duration taken at 986 kHz with PRFs of 200 Hz and 1000 Hz, respec- 
tively. The error bars in these and subsequent figures represent one 
standard deviation of the data above and one standard deviation below 
the mean threshold pressure. It should be noted that thresholds were 
undetectable at the 2 usec pulse width for the 200 Hz PRF, most likely 
due to lack of scattered signal from a two-cycle burst at such a low 
PRF. At 1000 Hz PRF, the thresholds are essentially constant down to 
the 10 usec pulse duration. At 5 and 2 usec pulse durations, the 
thresholds show a marked increase. The data in Figures 5.1 and 5.2 


show good agreement with Reference 6. 


Threshold Pressure (bar) 
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Threshold Pressure vs. Pulse Duration 986 kHz, 200 Hz PRF 
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Figure 5.2 
Threshold Pressure vs. Pulse Duration 986 kHz, 1000 Hz PRF 


Figures 5.3 and 5.4 are graphs of threshold pressure versus pulse 
duration at 2.17 MHz for PRFs of 200 and 1000 Hz, respectively. 
There is an essentially constant threshold at 200 Hz PRF for pulse 
durations of 15 usec and larger, while shorter pulse durations indicate 
a marked increase in threshold pressure. In contrast, however, at 
1000 Hz PRF, the thresholds are essentially constant for all pulse 
durations down to 2 usec, with a marked increase evident at a pulse 
duration of 1.5 usec. 

Figure 5.5 is a comparison of the thresholds obtained for the two 
source frequencies at 200 Hz PRF; Figure 5.6 is the same at 1000 Hz 
PRF. These results support theoretical predictions by Flynn [Ref. 4] 
and Apfel [Ref. 3] that the threshold is independent of frequency for 


sufficiently small cavitation nuclei. For parameters used in this 
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Figure 5.3 
Threshold Pressure vs. Pulse Duration 2.17 MHz, 200 Hz PRF 
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Figure 5.4 


Threshold Pressure vs. Pulse Duration 2.17 MHz, 1000 Hz PRF 
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Threshold Pressure (bar) 


Threshold Pressure (bar) 
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Figure 5.5 


Threshold Pressure vs. Pulse Duration 
986 kHz and 2.17 MHz at 200 Hz PRF 
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Figure 5.6 


Threshold Pressure vs. Pulse Duration 
986 kHz and 2.17 MHz at 1000 Hz PRF 
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experiment, threshold pressure values are essentially equal to Blake 
pressure thresholds, which have no frequency dependence (see 
equation 2.11). This conclusion is valud due to the close relationship 
between transient and Blake thresholds for the frequencies under 
consideration (see Figure 2.3). 

Figure 5.7 is a histogram of pressure thresholds (rounded to the 
nearest bar) measured for each individual transient cavitation event. 
Again from equation (2.11) above, there is a direct relationship 
between bubble size and Blake threshold pressure. Figures 5.8 
through 5.11 are histograms of threshold measurements for individual 
frequency/PRF combinations. It should be noted that the graphs may 
be slightly skewed due to the starting pressure amplitudes of most 
ramps. Since our experience indicated that expected threshold 
pressures were in the 6-9 bar range, source pressure ramps were 
started in the 3-4 bar range. An occasional low threshold event due to 
the presence of an unusually large nucleus may have been missed. 
Keeping this in mind, it appears that all histograms except Figure 5.8 
exhibit a roughly Gaussian distribution with mean values in the 7-8 bar 
range. This corresponds to bubble radii in the range of 0.075-0.09 
um. 

As stated above, when larger pulse width measurements are 
compared in Figures 5.5 and 5.6, the thresholds are independent of 
pulse duration. As pulse durations decrease, there is an increase in 


threshold. This observation would indicate that the larger pulse 
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Figure 5.7 


Histogram of All Threshold Pressure Measurements 
(rounded to nearest bar) 
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Histogram of All Threshold Pressure Measurements 
986 kHz, 200 Hz PRF 
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Figure 5.9 


Histogram of All Threshold Pressure Measurements 


986 kHz, 1000 Hz PRF 


20 





ara 
SAMs 


2 3 APS 6. re eerie 4 le ion) 7 12 > 70 2 227572425 
Threshold Pressure (bar) 


40 








30 








20 







Figure 5.10 


Histogram of All Threshold Pressure Measurements 


2.17 MHz, 200 Hz PRF 
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Figure 5.11 
Histogram of All Threshold Pressure Measurements 
2.17 MHz, 1000 Hz PRF 

durations allow time for a growth mechanism, such as rectified 
diffusion, to cause nuclei to evolve to a Blake or transient threshold 
radius. There appears to be insufficient time for the growth 
mechanism to take effect at shorter pulse durations. It is postulated 
that, during the rest time between pulses, nuclei dissolve to their 
original equilibrium sizes. 

Further examination of Figures 5.1 through 5.4 reveals that there 
is a distinct relationship between the flat and sloped portions of the 
graphs. Specifically, Figure 5.1 (low frequency/low PRF) reveals a 
continuously decreasing threshold pressure for increasing pulse 
duration. Although no additional data was taken at pulse durations 


greater than 40 usec, there were essentially no thresholds recorded at 
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levels lower than the threshold determined at that point, so it is 
anticipated that the graph flattens out beyond this point. An attempt 
was made to evaluate this hypothesis experimentally, but problems 
were encountered with the 986 kHz source transducer that prevented 
any further reliable data from being obtained at that frequency. Figure 
9.2 (low frequency/high PRF) shows relatively flat threshold 
measurements for pulse durations of 10 usec and greater, with 
significantly higher thresholds occurring at 5 and 2 usec pulse 
durations. Figure 5.3 (high frequency/low PRF) shows flat threshold 
measurements for pulse durations of 15-20 usec and greater, and 
increasing thresholds for 10, 5, and 2 usec pulse durations. Finally, 
Figure 5.4 (high frequency/high PRF) shows essentially flat threshold 
measurements down to 2 usec pulse duration. This observation 
prompted an additional set of 10 threshold measurements at 2.17 
MHz, 1000 Hz PRF, and 1.5 usec pulse duration (3 cycles per pulse) to 
ascertain whether there was a rise in threshold below 2 usec. The 
mean value of the ten readings was 12.9 bar, with standard deviation of 
2.90, which indicates a “knee” in this graph as well in the vicinity of 
the 2 usec pulse duration. There appears to be a relationship between 
frequency, pulse duration, and PRF which allows a certain range of 
pulses per unit time to mark the break point where shorter pulse 
durations begin to affect cavitation thresholds. In Reference 6 it was 
postulated that threshold became independent of pulse duration when 
pulses contained more than 10 cycles. Our measurements show the 


“knee” to be located as indicated in Table 5.1: 


L2 


TABLE 5.1 


PE N Vv. F ENCY 
PRF, AND PULSE DURATION COMBINATIONS 

f PRF Pulse length Cycles/pulse Cycles/sec 
(“KNEE”) PER ac 

Cycles/pulse) 
986 kHz 200 Hz 40 usec 40 8,000 
986 kHz 1000 Hz 10 psec 10 10,000 
2.17 MHz 200 Hz 15 psec 33 6,600 
Zag WillZ 1000 Hz 2 usec 4 4,000 


Although the value of PRF x (cycles/pulse) product varies from 
4,000 to 10,000, this analysis indicates that the knee of the graph is 
predictable within a reasonable range of pulse widths. In addition, our 
choice of pulse widths was not geared toward finding this break point 
of the graph. There may be other pulse durations that more accurately 
indicate the knee between those examined which yield a narrower 
range of cycles per unit time values. To test the hypothesis, a series of 
threshold measurements were taken at 2.17 MHz and 500 Hz PRF. It 
was anticipated that the knee of the curve should occur in the vicinity 
of the 10 psec pulse duration. Figure 5.12, which illustrates the 
results of these measurements, shows that the knee occurs as 
predicted. Although it is not prudent to generalize when discussing 
mechanisms as complex as bubble dynamics in sound fields, it is felt 


that this aspect of these findings bears further investigation. 
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Figure 5-12 
Threshold Pressure vs. Pulse Duration 2.17 MHz, 500 Hz PRF 


B. RECOMMENDATIONS 

All experimental measurements were taken using the passive 
acoustic scattering method, which proved to be an excellent detection 
method when pulse widths were large. The resulting signal scattered 
from the bubble “cloud” generated by the cavitation event resulted in a 
clear, distinct indication of a transient event. As pulse durations were 
decreased, however, the scattered signal was sometimes weak and 
lacking definition. It is essential that sonoluminescence detection 
and/or active acoustic detection schemes be fully incorporated into 
the procedures to minimize any systematic errors that may exist in 
the passive detection system. The optical detection system was 
limited only by unavailability of luminol (a sonoluminescence 


enhancer) until the very final weeks of the project. By then, it was too 
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late to bring the system to full operation. It is anticipated that this 
system, when fully operational, will allow the replacement of the 
human element of the current detection scheme with the electronics 
system. As gathering of cavitation data via manual observation of the 
transient signals on an oscilloscope requires exceptional vigilance and 
perseverance on the part of the operators, fully automating the 
detection mechanism is essential to the continued success of this 
apparatus. 

It is recommended, when the system is fully operational, that 
additional threshold measurements be taken at specific pulse dura- 
tion/PRF/frequency combinations in order to confirm the cycles per 
unit time effect noted in the hypothesis above. Extremely low or high 
PRFs (10 Hz, 8 kHz) should indicate whether the hypothesis has 
merit. Localization of a specific number of pulses or period of 
insonification which marks the knee of the threshold curve may lead 
to a better determination of the mechanisms which allow the 


cavitation events to fully evolve. 
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APPENDIX B 


F PHONE TRUCTION AND 
EN INF TION 


A. CONSTRUCTION DETAILS 

Our desire for small size led us to choose a large plastic syringe 
and a large-diameter hypodermic needle as the hydrophone housing. 
This gave us a convenient method for insulating the “hot” lead from 
the grounded elements, as well as providing ample storage area close 
to the active elements for the preamplifier circuitry. Our primary dif- 
ficulty in construction was working with small sizes of PVDF film. 
After numerous unsuccessful attempts to cut the PVDF to a shape 
which allowed for good electrical connection while minimizing active 
element size (1 mm2), we moved to a PVDF strip approximately 2 mm 
wide and about 4 cm in length. We reduced the bottom (positive) side 
of the film by etching away the PVDF using a Q-tip soaked in ferric 
chloride, which gave us a marginally accurate method of reducing 
element size while keeping a capability to have a reliable electrical 
connection. The top (ground) was not etched. The PVDF strip was 
then centered on a 3.6 mm x 6.4 mm brass plate attached to the “hot” 
lead. The barrel of the needle was insulated with several layers of heat 
shrink tubing to prevent any of the “hot” side from short-circuiting to 
ground. The remaining lengths of the PVDF strip were attached to the 
insulation along the needle barrel and connected to the surface of the 


needle with nickel paint. Gaps between the PVDF strip and needle 


ell 


were filled with epoxy adhesive to add stability. Resistance measure- 
ments were taken to ensure continuity of connections and that no 
short circuits existed. The unity gain preamplifier (voltage follower) 
was constructed and installed in the syringe and electrical connec- 
tions made to the hot lead and external ground connections. The 
entire apparatus was insulated by being coated with “Plasti-Dip.” See 


Figures B.1 and B.2. 


B. SENSITIVITY CALCULATION [Ref. 28] 


E = Vo/t where E = Electric Field developed 


Vo = Voltage Induced 
t = Thickness of Piezoelectric film 
Also, E = g33X3 where g33 is a constant indicating field/applied stress 


X3 = Applied mechanical stress 


Mo = Pressure Sensitivity = V/P = a a £33t 


Numerically, 


Mo = gaat = 339 x 10°3 —— (9 x 10-6 m) 


— A 
3.051 x 10-6 N/m = 3.051 x 10° 5 


20 logic (3.051 x 10-6) dB re 1V/Pa = -110.31 dB re 1V/Pa 


= 
5 
T 
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This calculation indicates that transducer pressure sensitivity is 
frequency independent. It should be mentioned that the calculated 
sensitivity is contingent upon the assumption that the nickel paint 
used to mount the PVDF to the brass mounting block effectively 
clamps the active face in the (1) and (2) directions, thus confining the 
face to operation solely in the thickness-expander mode. See Figure 


B.3. 
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Figure B.3 
PVDF Hydrophone Sensitivity Calibration Curves 


C. DESCRIPTION OF TEST AND CALIBRATION TECHNIQUE 

We chose the reciprocity technique for calibration of our 
hydrophone sensitivity because calibrated hydrophones in the 
megahertz range were not available for comparison. Our frequency 
band of interest was 0.5 MHz to 4.0 MHz. For the reciprocity 


calibration, we used two Panametric V323 transducers as the trans- 
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mitting and reversible elements of the experimental setup. These 
transducers are designed for operation at 2.25 MHz and provided 
smooth response through the frequency range. The calibration set up 
was as noted in the equipment diagram included as Figure 4.1. The 
distance between the transducers was 30 cm, which is well into the 
far field for all frequencies of interest. Far-field calculation for the 


maximum frequency is as follows: 
L2 
r=77 [Ref. 25] 


where L = 5 x 10-8, m = largest extent of source active element, A = 
wavelength at 4.0 MHz in fresh water, r = minimum extent of far field, 


and 


ig ee = 0.1689 m = 16.89 cm 
We took Voltage Level (in dBV) and input current (in milli- 
amperes) readings from the HP 3585A Spectrum Analyzer and HP 
3478A Digital Multimeter in accordance with standard reciprocity 
techniques. 
Open-circuit voltage sensitivity of our PVDF hydrophone was 


determined using the following equation: 


oN 2 
De VxXVx 
Mo, = iq WaIE (Ref. 25] 


When converted to sensitivity level, this calculation became: 
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| 2 
ML = Vig a ~Vr,— IL, — 20 log ime 20 log 2 x5 


D. DIRECTIVITY EFFECTS/EXCITATION OF FIRST THICKNESS 
MODE 


1. Analysis of Directivity Effects 

One possible cause of the difference between actual and pre- 
dicted sensitivities may be the change in the directivity of the 
hydrophone as frequency increases. As the transmitting transducer 
frequency is increased, both the hydrophone and transducer become 
more directional, implying more efficient transmission by the trans- 
ducer and more efficient reception by the hydrophone of the acoustic 
power in a given direction. The trade-off is that the main lobes of 
both the transducer and the hydrophone become narrower. Using the 


approximation for a piston-like source [Ref. 25], 


D = 1/4 k2L)L92 where L) = 2 mm and Ly = 5 mm 
fe= 7c = 2nt/C 
At 500 kHz, DI = 10 log D = 10.5 dB 


These “corrected” calculated values are plotted as such in Figure B.3. 


2. Iculation of First Thickn M of Brass Mountin 
Plate 
It was originally believed that the brass mounting plate would 
behave as if it were perfectly rigid in the frequency range of interest. 


Measurements led us to believe that this may not be the case. For this 


Su 


reason, it was decided to determine if the first thickness mode of the 
brass mounting plate was being excited somewhere in the range 
between 0.5 MHz and 4.0 MHz. The first thickness mode will be 
excited with t = 4/2 or when A = 2t. For a thickness of 1.3 mm, A = 
2° mm. Smice c = it =e 


f=c/A = (4700 m/s)/2.6 x 10-3 m = 1.81 MHz 


It becomes evident that it is likely that the first thickness 


mode is being excited to some extent. 


E. DISCUSSION 

As noted above, theoretical sensitivity was calculated assuming 
that the longitudinal and transverse modes were clamped and that 
only the thickness mode was excited. Data showed little agreement 
with these original theoretical calculations, however. The PVDF 
hydrophone’s active face was considerably larger than the wavelengths 
of the frequencies of interest, which makes the receiving beam pat- 
tern more narrow as the frequency increases. This directivity also 
adds considerably to the sensitivity when sound is received “on axis.” 
Directivity and Directivity Index as calculated previously in this paper 
contributed greatly to the sensitivity of the hydrophone. As can be 
seen from the graph of Pressure Sensitivity vs. Frequency (See Figure 
B.3), there is good agreement between theoretical and experimental 
values at the lower end of the frequency range (up to 1.5 MHz). At 


higher frequencies, additional variables are introduced which limit 
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agreement. A likely contributor to the higher frequency disparity 
might be slight misalignment of the hydrophone. As frequency goes 
up, the beam pattern becomes increasingly narrow. Slight misalign- 
ment could result in poorer reception of the transmitted acoustic 
energy. Lastly, as alluded to previously, it is possible that the first 
thickness mode of the brass mounting plate is being excited in the 


frequency range of interest. 
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